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Abstract 

We study the spectrum of off-diagonal fluctuations between displaced fuzzy spheres in 
the BMN plane wave matrix model. The displacement is along the plane of the fuzzy 
spheres. We find that when two fuzzy spheres intersect at angles classical tachyons 
develop and that the spectrum of these modes can be computed analytically. These 
tachyons can be related to the familiar Nielsen-Olesen instabilities in Yang-Mills theory 
on a constant magnetic background. Many features of the problem become more 
apparent when we compare with maximally supersymmetric Yang-Mills on a sphere, 
of which this system is a truncation. We also set up a simple oscillatory trajectory on 
the displacement between the fuzzy spheres and study the dynamics of the modes as 
they become tachyonic for part of the oscillations. We speculate on their role regarding 
the possible thermalization of the system. 
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1 Introduction 

In gravitational theories it is expected that very high energy scattering at small impact 
parameter is dominated by the formation and evaporation of black holes. Our basic under- 
standing of how this might happen in a semiclassical approximation has led to Hawking's 
paradox, the statement that black holes might lose information with the scattering being 
non-unitary [1] . The same understanding that has led to this paradox shows that black holes 
have a very large entropy and that they also have a temperature. Both of these features 
cannot be explained by classical means. 

The gauge/gravity duality has given us a way, in principle, to formulate these phenomena 
in a unitary framework, where the black hole intermediate state should be described by some 
approximately thermal object with large entropy in a dual quantum field theory. The details 
of such approximately thermal field theory configuration, however, have not been worked out 
so far. Large black holes in AdS have positive specific heat [2] and can be readily identified 
with a thermal state in the field theory [3]. On the other hand, small black holes are not 
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stable objects (since they evaporate) and a complete description of the system would require 
a theory of how such small black holes form and evaporate. An initial state that produces 
the black hole is not thermal, but the dynamics should be such that the system thermalizes 
rapidly (the black hole formation and ringing suggests this kind of picture The ensuing 
evaporation of the black hole shows that the system should be thought of as being out of 
equilibrium, but very long lived. For this second stage of evaporation one would need a 
detailed analysis of how the degrees of freedom - the entropy - of the black hole escape from 
it, if one is to claim to have really solved the black hole information paradox. 

An attempt to describe the initial stages of thermalization of small black hole^j]was done 
for M = 4 super Yang-Mills theory in [6], but the details of the dynamics in that context 
are still poorly understood. In general, problems of thermalization in full-fledged quantum 
field theories are hard to approach, especially at strong coupling. A promising strategy to 
try to address these issues is therefore to look for examples with fewer degrees of freedom. 
This paper is a first step in this direction, where we show that formulating the problem of 
black hole thermalization in a simpler setting has various advantages]^] One of such settings 
is represented by matrix quantum mechanical systems. 

The BFSS matrix model j8] is the prime example of a matrix quantum mechanics. It is 
dual to M-theory on a discrete light-cone quantization of flat space and it has been argued 
that it accommodates black holes, which have been analyzed in various works [9]. The 
geometry of the black hole horizon can be related to properties of the matrix quantum 
mechanics and to the presence of tachyons [ID] . Some of the studies of the BFSS model 
have also included numerical analysis based on lattice methods in Euclidean field theory 
However, a detailed analysis of the formation and thermalization properties of these 
matrix black holes is still missing. Moreover, a numerical instability arises in these Euclidean 
thermal systems, exactly because eigenvalues can leave the black hole state. As such, the 
Euclidean ensemble does not exist, as one would expect when studying systems that cannot 
be in true equilibrium. 

This problem can be fixed by considering instead the BMN matrix model [12], where 
the large volume instability of having a moduli space of vacua with a continuous spectrum 
is cured and the simulations make sense [13J. The BMN matrix model is in fact a mass 
deformation of the BFSS model which lifts the flat directions and has no problems from this 
point of view. As the BFSS model, it has an M-theory interpretation, as it describes the 
discrete light-cone quantization of the theory on a plane wave (rather than on flat space). 
Another nice feature of the BMN matrix model is that it can be thought of as an SU(2) 

1 These are black holes which are much smaller than the radius of the AdS geometry and behave as 
ten-dimensional Schwarzschild black holes [S]. 

2 Other studies of thermalization based on toy models have been carried out in [TJ, but part of the system 
there already begins in a thermal state and the process of black hole production is not addressed. 
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invariant sector of M = 4 super Yang-Mills theory on S 3 x 1R p3] . This means that certain 
classes of analysis of this model might also be fruitful to better understand the physics of 
that four-dimensional theory. 

Our main interest in this paper is to explore some time dependent dynamics, in a matrix 
quantum mechanics setting, that might help us understand how the initial stages of black 
hole thermalization might proceed. For this, we need some simple configurations with given 
initial conditions that are similar to the problem of scattering gravitons at high energies. We 
can then study more precisely the dynamics of the modes that produce the thermalization. 
In the analysis of [6] , it was argued that particle production in off-diagonal modes connecting 
eigenvalues was responsible for generating the entropy of the black hole configuration and for 
trapping the eigenvalues. If we try to mimic this intuition, we would expect it to work very 
similarly for the BFSS matrix model: we could think of scattering two graviton states with 
large matrices and to hopefully see the black hole form. Unfortunately, we are unable to do 
this because these graviton states in the BFSS matrix model are very poorly understood: 
they are bound states at threshold and therefore we would need a detailed understanding of 
their wave function, which we do not have, to do the analysis. 

For this reason and for the reasons mentioned above, we turn instead our attention to 
the BMN matrix model. Here the set of ground states is richer than in the BFSS model and 
instead of scattering eigenvalues off each other we can scatter more complicated ground state 
configurations. These configurations are given by sets of concentric fuzzy spheres. Each such 
fuzzy sphere is a giant graviton [15] that can grow in size due to the Myers effect |16j . For 
each fuzzy sphere, one has a decoupled U(l) set of degrees of freedom that correspond to the 
center of mass motion of the corresponding membrane object. These degrees of freedom can 
be excited leaving the rest of the SU(N) degrees of freedom unaffected. This means that the 
fuzzy spheres can be rigidly moved around without being deformed and their geometry is 
simple to analyze. Thus, it is a simple matter to find reasonable initial conditions for these 
fuzzy spheres. By aiming them properly, we can make them collide (namely, have two such 
fuzzy spheres intersect each other at some time t). The dynamics near the intersection is 
interesting. From the point of view of string theory, these are configurations with branes at 
angles, so there is a possibility of having tachyons form at these intersections [T7] . Moreover, 
even in the absence of tachyons, modes near the intersections are light. As we move the 
branes past each other, these modes have a time dependent mass and in general we expect 
copious particle production in these modes if the adiabatic approximation breaks down. At 
weak coupling this is similar to the problem of preheating in inflation [18]. In our case, 
the matrix model black hole would result from thermalization after moduli trapping [T9]|^] 
A further advantage that we have in this setup over the BFSS model is that the classical 
solutions corresponding to branes crossing each other are periodic: the branes can cross 

3 The trapping between DO-branes was analyzed in |20j . 
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repeatedly and the fluctuations that are produced in this way have time to grow through 
many repetitions until they cause a large back-reaction. 

The main goal of this paper is to describe in detail the mass spectrum of the bosonic 
modes stretching between fuzzy spheres, as we displace the spheres along a direction longitu- 
dinal to their world- volume and we make them intersect. Our calculation is a generalization 
of the analysis done in [2TJ to study the fluctuations around concentric fuzzy spheres]^] While 
fluctuations around some fuzzy sphere configurations have already been considered in the 
past by other groups, in this paper we focus our attention on a different and novel set of con- 
figurations. In previous studies, the displacements of the fuzzy spheres are in the transverse 
directions to the spheres and the distance between the spheres (or other branes [21]) does 
not change in time, or does so in such a way that the time dependence of the off-diagonal 
modes is perturbatively small [25J. 

We will argue that in our setup there is a spectrum of tachyons at the intersection 
locus of the fuzzy spheres. The intersection locus being one-dimensional gives us a one- 
dimensional tower of such tachyons. There is such a tachyon even for two DO-branes at some 
finite distance. These tachyons can have large negative mass when at least one of the fuzzy 
spheres corresponds to a large matrix, so that any small fluctuation in these modes can grow 
rapidly and might lead to fast thermalization. A detailed analysis of the evolution of the 
system following the tachyon production is beyond the scope of the present work, but we 
sketch nonetheless a picture of what happens in a simple case. 

The paper is organized as follows. In Section [2] we review briefly the BMN plane wave 
matrix model and the classical vacua of the theory. We present a novel derivation of the 
spectrum of fluctuations around concentric fuzzy spheres making use of the aforementioned 
relation between the BMN matrix model and M = 4 super Yang-Mills. In Section [3] we 
consider a classical solution of the BMN matrix model consisting of two fuzzy sphere which 
are not concentric, but displaced along a direction longitudinal to their world-volume. After 
giving a geometrical characterization of the fluctuations around this configuration, we pro- 
ceed in Section [4] with a detailed computation of the spectrum of off-diagonal fluctuations, 
both along the transverse directions to the spheres and along the longitudinal ones. We 
find that certain longitudinal modes can become tachyonic. In Section [5] we study the time 
dependence of these tachyonic modes as we allow the displaced fuzzy spheres to oscillate 
one towards the other along the direction of the displacement. Some concluding remarks 
can be found in Section [6j while in the Appendix we present an alternative derivation of the 
spectrum of longitudinal fluctuations. 

4 Some studies of the thermodynamics of these configurations include the works in (22]. Other important 
result about protected multiplets can be found in [23j . 
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2 The BMN matrix model 



The BMN matrix model [12] is a massive deformation of the BFSS matrix model [8]. The 
latter is obtained from the dimensional reduction of ten-dimensional M = 1 super Yang-Mills 
down to + 1 dimensions and has an action given by 



Sbfss = ^ J dt ti (^D t X r ) 2 + ^[X 1 , X J ] 2 ^j + fermions , 



(2.1) 



where X 1 (I = 1, . . . , 9) are nine hermitian matrices. The covariant time derivative is given 
by 

D t X T = d t X* -i^X 1 ] (2.2) 

and g is a dimensionful coupling constant that can be removed by rescaling the fields and the 
time coordinate. It can be set to one, if desired, or factored out of the action and interpreted 
as determining h. We will not work in detail with the fermions in this paper, so we shall 
just suppress them from now on. The Hamiltonian of this system (in the A t — gauge) is 
given by 

H= 1 -tr^r-l- 2 [X I ,X J ]^ . (2.3) 



The BMN matrix model system is a massive deformation of (2.1) that preserves all 32 
supersymmetries. It also preserves a diagonal set of modes that decouple and constitute a 
system of free degrees of freedom. These are the 'center of mass motion' degrees of freedom in 
the BFSS matrix model. The BMN matrix model splits the X 1 into two groups of variables: 
X 1 ' 2 ' 3 , that we will label X 1 , and X 4 ' " ' 9 , that we will label Y a . The action includes additional 
terms given by 

Sbmn = Sbfss - ^ J dttr (V(X 1 ) 2 + ^(Y a ) 2 + 2 / j,ie ijk X i X j X k ^ . (2.4) 

In the conventions above, \x is real and has been rescaled by a factor of 3 with respect to 
[12] . It has units of frequency, as X 1 . The equations of motion following from this action 
are 

X 1 = -n 2 X i -3i t i^ jk X j X k -[[X i ,X I ],X I ] , 

Y a = ~^Y a - [[Y a ,X J ] ,X J ] . (2.5) 

It is convenient for our study to recast the potential for the X % fields in the following 
form 



V (X) - — tr 

v BMN — 2g2 LI 
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{i[X 2 ,X^+ l xX x ) 2 + {i{X\X 1 ]+ l xX 2 ) 2 +{i{X\X 2 ] + ^) 2 . (2.6) 
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In this paper we will choose to rescale — > 1 by both scaling the matrices and the time 
coordinate. This can always be done. The overall action then has a factor of 1/g 2 in front. 
This serves as a calibration of h: we are free to absorb g 2 in the definition of h. This has 
no effect on the classical physics except for the global normalization of the energy units. 
It is only when we quantize the theory that the value of g will be important and it will 
characterize the strength of quantum effects or fluctuations around classical solutions of the 
dynamics. 

An important point to notice is that the potential is a sum of squares of hermitian 

(X) 

matrices. The VgjJ N is positive definite and the same is true for the other terms in the 
potential: the quadratic terms in Y a are obviously a sum of squares and the rest of the 
terms in (2.1) are commutator squared terms with the right sign to make them positive 
definite. 

Another very useful result to recall is that the BMN matrix model can be obtained by 
considering a truncation to the 577(2) invariant configurations of M = 4 super Yang-Mills 
on S 3 x R [U\. The sphere has a 50(4) ~ SU(2) L x SU(2) R symmetry and the SU(2) L 
invariant sector of the theory gives exactly the model above. In that case, g is identified 
with the coupling constant of A^ = 4 super Yang-Mills. The X % degrees of freedom arise 
from the gauge connection on S 3 and in the usual field theory analysis would have mass 2. 
The Y a degrees of freedom arise instead from the scalars <p a of the super Yang-Mills theory 
and would ordinarily have mass equal to one. This corresponds to setting fi = 2 in the BMN 
model above. Equivalently, we can think of this as quantizing the M = 4 theory on a sphere 
of radius equal to 2, rather than one. This SU(2)l invariant reduction is a convenient device 
to calculate sometimes and it also explains why the potential is a sum of squares. 



2.1 Fuzzy sphere ground states and fluctuations 

The ground states of the BMN matrix model are those that have zero energy. These are 
characterized by Y a = and also by 

[X i ,X*]=ie ijk X k , (2.7) 

where we have set \x = 1. These are obtained by requiring that each of the individual squares 
in the potential vanish. 

The most general solution to this equation is given by a possibly reducible lie algebra 
representation of SU(2), where X % ~ © a -^(n )> with the sum indicating a sum over irre- 
ducible representations of SU(2). Each irreducible representation (which we can label by 
the size of the representation n = 2j + 1, or by the maximal spin j) gives rise to a fuzzy 
sphere configuration. In the BMN model these are interpreted as giant graviton membranes 
of the plane wave limit M-theory dual. 
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If we fix the total size of the matrices N, we can set N = ^2 a n a . The complete set 
of vacua of the theory is characterized by these possible splittings. This is equal to the 
partitions of N. 

Given such a configuration, we can think of the ground state vacuum expectation values 
as being block diagonal. We can then ask what is the spectrum of off-diagonal fluctuations 
connecting two of these fuzzy spheres. This has been analyzed in detail in [21 J . 

Here, we will reproduce that result using a slightly different calculation. The idea is to 
remember that this matrix quantum mechanics is an SU(2) l invariant reduction of the U(N) 
M = 4 super Yang-Mills theory on S 3 x K. All of the fuzzy sphere vacua are related to each 
other by a gauge transformation in the theory on S 3 [26]. The gauge transformation that 
relates them is not SU(2)l invariant, so the modes that are considered SU(2)l invariant 
get shuffled. The gauge transformation uses the obvious map S 3 — » SU(2), since they are 
identical spaces. This SU (2) has to be embedded in the gauge group. 

For a fuzzy sphere configuration N = X]a n a> we embed SU(2) into U(N) by the action 
on the fundamental of U(N) according to © a -R a = ® a (n a ), where we are labeling the 
representations by U(N). This SU{2) embedding can also be thought of as SU(2) L : the 
translation on the sphere by SU (2) l generates the sphere itself. Also notice that the fuzzy 
sphere configurations are SU(2)i invariant only if an SU(2)l rotation is accompanied by a 
compensating gauge transformation. 

Now let us do some group theory. The Fock space spectrum of physical polarizations of 
fluctuations of the A = vacuum is given by 

Spec(4>) = 

j 

Spec(A) = @(j,j + l)©(j + l,j), (2.8) 
j 

where we are using the spin notation for the representations. The first line indicates the 
spectrum of representations of the scalar fluctuations, while the one of the second line rep- 
resents the transverse fluctuations of gluons. The energy for the scalars (j) is 2j + 1, while 
for the vectors it is 2j + 2. 

When we consider the embedding on matrices, we have that the off diagonal block con- 
necting block m and 1112 transforms as j\ ® 32 with respect to SU(2). This is, the matrices 
transform as Of^^-j^f with respect to SU(2)l- If we tensor these together, we find that 
for the scalars we need to take the tensor product 

0./ (./•./) -0 U':l) ■ (2-9) 

j j j'=\j-j\ 
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We can only have SU(2) L singlets if / = 0, so we find that the singlet sector is given by the 
case where j = j. We find this way that 

j j 

with energy 2j + 1. Notice that j is either only half integer or only integer, as per the usual 
rules of angular momentum. 

Similarly, we find that for the vectors 

j ® + 1) © (j + 1J)) \ singlet * 0(0, j + 1) © (0, j - 1) , (2.11) 

j j 

with energies 2j + 2 and 2j respectively. For the special case j = 0, the representation 
(0,-1) is not counted. 

Properly normalizing to the value of \x we have used (i.e., dividing by 2 the energies 

above), we get that on the fuzzy sphere the scalars Y a have a spectrum under SU(2)r 
rotations given by 

%c(y) = 0], (2.12) 
3 

with energies j + 1/2, where j is only integer or half integer and it runs between \ji — J2I 
and ji + J2 in integer steps. 

Similarly, for the X 1 fields, we get a group decomposition 

Spec(X) = 0(j + 1) © (j - 1) , (2-13) 

3 

with energies j + 1 and j respectively. Remember that j is always integer or half integer. 
The same analysis can be done for fermions, we will not do this here. 

The representations j appearing in the decomposition of the off-diagonal matrix fluctua- 
tions of Y a are called fuzzy monopole harmonics. They are fuzzy spherical harmonics if the 
two representations have the same dimension. 

For the X 1 variables, the j representations appearing in the decomposition are called 
fuzzy monopole vector (tensor) harmonics and fuzzy vector (tensor) harmonics if the two 
representations have the same size. 

Notice that we only described the physical fluctuations. The typical configurations also 
have zero modes due to gauge transformations. These zero modes have been projected out. 
If we add them, we get that X has an additional set of zero modes that are described by the 
fuzzy monopole spherical harmonics. There is one more thing that needs to be remarked: 
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the zero mode associated to rotations in the diagonal U(l) is absent as nothing is charged 
under it. 

For matrices of the same size the theory has an enhanced SU(2) unbroken gauge symme- 
try and the modes organize themselves into triplets and singlets of SU(2). This generalizes 
to when one has more coincident fuzzy spheres. 



3 Kicking the spheres 



We now consider the following classical solutions of the BMN matrix mo delfl 

...\ 



/L\ ni) + Re(b\*l {ni) eMit)) 







> 2 ) +^e(6 J 2 * l (na) exp(it)) ... 



V 



(3.1) 



We have turned on modes proportional to the identity that decouple on each diagonal block 
and describe the center of mass motion of the fuzzy spheres. These modes leave the shape 
of the spheres invariant, but move their centers of mass in time. A similar kicking of the 
spheres can be done by turning on a displacement b a a along the Y a directions 



on 



fUe(bf * l (ni) exp(it/2)) 

3te(&g * l (n2) exp()'t/2)) 

V 



• A 



(3.2) 



Notice that in order to obey the equations of motion (2.5) the Y a directions must have a 



frequency that is half the frequency of the X 1 directions. For the rest of this paper we shall 
set (Y a ) = 0. 

We restrict now our analysis to the case of two diagonal blocks, with ranks ri\ and n%. 
We will first study the dynamics of the off-diagonal modes connecting the two fuzzy spheres 
when they are large, i.e. with m, n-2, and n% — n 2 being large. For the configurations with 
b\ = 0, the off-diagonal modes have a high angular frequency starting at n\ — n 2 , so this 
frequency is much larger than the frequency of the zero mode that we are kicking. This 
means that small fluctuations on these degrees of freedom are expected to be adiabatic for 
a typical b l a . 

There is a clear procedure to analyze this type of configurations. One thinks of this 
as a Born-Oppenheimer approximation where one can solve the dynamics of fast degrees 
of freedom as we freeze the slow degrees of freedom. The fast degrees of freedom are the 
off-diagonal modes and the slow degrees of freedom are going to be the zero modes that we 



5 These configurations are non-BPS. 
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turned on. They are only fast relative to the motion we have described if their spin is large. 
We will analyze these by freezing the result at some time t and studying the spectrum of 
quadratic fluctuations of the off-diagonal modes for that frozen configuration. Since t is fixed 
for the analysis, we can use the rotational symmetry of the system to align b\ — 62 along the 
3-axis {i.e. only b\ — b\ 7^ 0). Moreover, a combination of the form ni5i+n2&2 is proportional 
to the center of mass of the whole system and is a decoupled degree of freedom. This means 
that all the dynamics we are interested in depends only on b\ — b 2 . We can without loss of 
generality set b± = 0. Also, we can choose the displacements to be real. 

In the rest of this section we will give a geometric characterization of what type of results 
we expect to find. A detailed analysis of the spectrum of these fluctuations will then follow 
in the next section. 

3.1 Geometric characterization of the dynamics 

These matrix configurations can be described geometrically in the M-theory plane wave ge- 
ometry. A fuzzy sphere of rank n corresponds to an M2-brane giant graviton of size propor- 
tional to its light-cone momentum n (the proportionality constant depends on conventions). 
In matrix coordinate units, n acts as a cut-off on the range of the angular momentum (see 
|27j). so that the size of the sphere is also proportional to the maximum eigenvalue j. Defin- 
ing the radius of the sphere as given by the distance to the center of mass (this is also done 
in Matrix Theory [28J), we get 

R 2 = (X 1 - x\ m ) 2 + (X 2 - x 2 cm ) 2 + (X 3 - x 3 cm ) 2 = 3(3 + 1) . (3.3) 

At large j we have that the radius is R = j + 1/2 + o(l/j). Notice that the spacing of the 
radii between different j is essentially Aj. These M2-branes are described as D2-branes. Off- 
diagonal modes connecting different such configurations are interpreted as strings stretching 
between the D2-branes, while the eigenvalues are interpreted as DO-branes. It is standard to 
think of the D2-branes as branes that have absorbed DO-branes and as a consequence have 
a strong magnetic field on their world- volume [29|. As a matrix of rank n has n eigenvalues, 
this is the magnetic flux threading the D2-brane sphere. A string ending on the D2-brane 
is charged under this magnetic field and experiences a magnetic monopole flux of strength 
n. As is well known, if we consider a charged scalar degree of freedom on a sphere with 
a magnetic monopole background and if we restrict to the lowest Landau level, the wave 
functions carry angular momentum and can be argued to be localized on the sphere so that 
they cover an area of order 1/n. The angular momentum points in the direction on which we 
localize the wave function on the sphere. If we reverse the charge, the angular momentum 
points in the opposite direction. 

If we have two fuzzy spheres of sizes j,j' at rest they can be described by two concentric 
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spheres of ranks n,n'. Now let us consider strings stretching between them. If we attach n 
states to one sphere and n! states to the other sphere, we get a total of nn' possible strings. 
These strings will have an endpoint on a sphere that is associated to a positive charge and the 
other end on the other sphere will have the opposite charge (the string theory is oriented). 
The effective magnetic flux that the particle sees is n — n' . This describes the minimal 
angular momentum that the modes connecting the fuzzy spheres can have: \j — This is 
also the length of a string stretched from the north pole of one fuzzy sphere to the north 
pole of the other one. 




Figure 1: An off-diagonal mode can be thought of as a string stretching between two D2- 
branes. In this figure we consider the case of zero displacement, i.e. concentric spheres. 
The angular momentum of a state is given by subtracting the position vectors of the string 
endpoints. Using the symmetries of the system, the angular momentum vector can always 
be chosen to be aligned along the vertical axis. 

The string of maximum length has angular momentum given by j + f and a mass of 
order j + j' in appropriate units. The angular momentum vector can be obtained by taking 
the difference of the position vectors of the endpoints of the string. This is depicted in figure 
[TJ The angular momentum vector of the string state points parallel to the string. Also, the 
geometric length is proportional to the length of the angular momentum vector (this also 
happens for Non Commutative field theories on the Moyal Plane [30]). We depict in the 
figure various highest weight states. The longest string goes from the north pole of one fuzzy 
sphere to the south pole of the other. Notice also that each string endpoint can be thought 
of as occupying a uniform fixed area on each sphere. The sphere with n\ eigenvalues has 
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rii such patches and similarly the second one, with eigenvalues, will be divided into ri2 
patches. Each of these is to be though of as a DO-brane end-point (region) on the sphere. 

If we compare with the spectrum of the Y a fluctuations, we get a precise matching 
between the possible values of angular momentum we compute geometrically in this way and 
the ones we obtain from the field theory calculation. These are transverse polarizations of 
the strings to the three directions in which the branes are embedded. For polarizations of the 
string modes in the brane 3-plane, they have extra spin: indeed, they carry one unit of spin 
that is either along the direction of the string, or opposite to it (only transverse polarizations 
appear on the string) and one can match this to the values of angular momentum of the X 1 
fluctuations as well. Again, the geometric estimate of the mass is good enough. Notice that 
there are in general correction of order one to the mass. 

Now we can consider what happens when we displace the spheres (see figure [2|. It is 
clear that the density of the string endpoints on each fuzzy sphere is not going to change. 
This is because this is roughly the density of eigenvalues per unit area. Moreover, since the 
end points of the string are charges in a very strong magnetic field, the magnetic field is 
good at keeping them from changing positions. Roughly, we should imagine that the string 
ends are stuck to a particular DO brane and that these DO branes are heavy and don't get 
moved around much by the force of a string. Indeed, in perturbation theory the tension 
of the D-branes is large [31]. So we can argue that the endpoint location of the string on 
each fuzzy sphere will not change, nor how we think about its angular momentum in the 
3-direction (the one that is preserved by the configuration). However, the length of the string 
will change. 




Figure 2: The length of the strings changes as we displace the fuzzy spheres. 
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If we use the same labels for the string endpoints as before, that is we label them by 
their angular momentum, we find that the length is given by 

L 2 ~ (AL 1 ) 2 + (AL 2 ) 2 + (AL 3 - bf , (3.4) 

where b is the displacement. The masses should then be roughly given by 

M 2 ~ (AL) 2 - 26AL 3 + b 2 . (3.5) 

The mass formula will attain the minimum value on a sphere for fixed (AL) 2 when AL 3 
takes either the maximum or the minimum value, depending on the sign of b. 

Clearly, this value will be minimized when the spheres touch. The states with minimum 
energy have their spin aligned along the 3-axis shown in the figure (the strings of length 
zero do not have an horizontal component of L). There will be corrections to this simple 
geometric formula. This can be seen from comparing the values of the energy when the 
spheres are concentric to those that are obtained from the exact answer. These corrections 
to the mass squared are of order j and can in principle make some of the modes tachyonic. 
This is what one can expect from the intuition of branes at angles [17J. To check this requires 
doing a computation, that we do in the next section. What is important to notice is that if 
there are tachyons, they are localized near where the spheres touch. In this examples this 
is a circle, so one can expect that these degrees of freedom are like a one-dimensional field 
theory set of modes. 

Notice that the naive geometrical correction is proportional to AL 3 . If we treat it as 
an operator in the theory of angular momentum it commutes with (AL) 2 , so even though 
the configurations break the rotational symmetry, the spherical harmonic representation of 
the states should remain diagonal. This is similar to what happens in the computation of 
the Zeeman splitting in the hydrogen atom. We will see this explicitly when we do the full 
computation. Also notice that the most tachyonic mode will depend on the displacement, 
because of all the vertical strings the one that is the shortest depends on the precise value 
of b. This means that since in our kicked sphere solution the displacement is changing with 
time, which off-diagonal mode is the most tachyonic depends on time. If all of these start 
condensing, one expects that, since they carry angular momentum, the axial symmetry will 
be broken classically: the spheres will deform non-uniformly and change shape. The L 3 
will remain a constant of motion and the axial symmetry will only be restored quantum 
mechanically by averaging over the orientation of the resulting shape. The axial symmetry 
can be restored later again classically if the system thermalizes and any coarse grained 
observation becomes sufficiently homogeneous (this usually requires large N). Such effects 
would be expected when we form a non-rotating black hole in the gravity theory. 
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4 Spectrum of fluctuations 



Now we are ready to start calculating the spectrum of fluctuations of off-diagonal modes 
for displaced fuzzy spheres. As we have described above, at large N the off-diagonal modes 
connecting two fuzzy spheres are generically heavy, except where the fuzzy spheres intersect. 
We have also argued that geometrically it seems that the basis of spherical harmonics is 
preserved. 

We will expand both the X 1 and Y a variables in spherical monopole harmonics. This is 
because the off-diagonal matrices have the same rank. The fuzzy vector spherical harmonics 
will be particular linear combinations of these. The basis of matrices will be labeled by 
Ytm £ Hom(n2,ni), which explains how we think of these as matrices. The conjugate 
harmonic is Y} m = (— l) m Yt_ m e Hom(ni,n 2 ). They are normalized so that 

tr(F/ m Y£/ m /) = ^6 U 'Smm' ■ (4.1) 

The matrices Y a that are off-diagonal and connecting the two fuzzy spheres will be 
hermitian and expanded as follows 

Y a = vLYim + (VLTYL , (4.2) 

e,m 

or, writing the blocks more explicitly, 

t 




m)*^m 



(4.3) 

We will do the same for the X 1 matrices. The difference is that the X % matrices will also 



have a vacuum expectation value as in (]3.1 ). As argued above, we can take 

o 

3 



<* 3 > = ( <"> ; |, (4.4) 



where b is the displacement of the second fuzzy sphere, while we take (X 1 ) and (X 2 ) with no 
displacement. Our goal is to compute the masses of the fluctuations in terms of £, m, ni, n 2 , b 
and any possible mixings between the states. To do this, we will need the following commu- 
tation relations for the Y^ m 



[L 3 ,Y lm ] = mY t 



[L + ,Y em ] = y/(£-m)(e + m + l)Y em+1 , 

[L-,Y em ] = v /(£ + m)(4? + l-m)y £m _ 1 , (4.5) 
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and their adjoints 



[L + ,Yl\ = - v /(£ + m)(£-m + l)r/ m _ 1 , 

[L-,Yt\ = - v /(£-m)(£ + m + l)r/ m+1 , (4.6) 

where = L 1 ± iL 2 . 

First we study the fluctuations of Y a , which are easier to analyze, and then those of 
X % . We will do the calculation for arbitrary values of ni,n2 with the configurations being 
off-shell (frozen at finite displacement). 



4.1 Transverse fluctuations 

To compute the spectrum of transverse fluctuations we need to expand the Lagrangian in 
fluctuations 5Y a = J2 e m Sy^ m Yi m + (Sy^ m )*Y^ m keeping only the terms up to quadratic order. 
Since the Y a do not have a vev and they appear always at least quadratically in the action, 
we can just truncate the action to quadratic order in Y a and replace the X 1 by their vevs. 
There can be no mixing between X 1 and Y a (this is easiest to argue by the fact that the Y a 
carry SO(6) symmetry labels and the configurations under study are 5*0(6) invariant). The 
expansion of the Lagrangian is straightforward and yields 

£ (y) = hr(Y a ) 2 -hr(Y a ) 2 + hi[(X l ),Y a ] 2 
2 o 2 

= \^(Y a ) 2 - \tr(Y a ) 2 + ^tr[(X 3 ), Y a ] 2 + 1 -tx[(X+), Y«] [(X~), Y a ] . (4.7) 
The first thing we compute is the kinetic term, which is obviously given by 

42 = ^Ei^i 2 - ( 4 - 8 ) 

£,m 

For b = it is easy to see that the potential can be written as [21] 



£ (y) = --tr 

^mass 2 



ya [L\ [L\Y a ]]^j 



(4.9) 



Since [L l , [L\ Ye m ]] = t(t + l)Y^ m , this implies that at b = we get that the mass terms are 
given by 



This matches of course with the result obtained in section 2.1 using group theoretical argu- 



ments. A more detailed description of these spherical harmonics can be found in 
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We now turn on b. Notice that schematically it is (X 3 ) = L 3 + b 




1 



so that 



[{X 3 ),5Y a } = [L 3 ,5Y a }+b 




,8Y a 



IJm - b)5y a lm Y lm - (m - b){5y1 m YY) 



tin 



(4.11) 



Notice that the addition of b does not mix the Yi m with each other in the commutator and 
all it does is to replace m — > m — b in the commutation relations. This is what makes 
the computation so simple. This is what we were arguing for geometrically in the previous 
section. When we square the expression above, we get the same result as when 6 = 0, plus 
additional terms that are cross terms and a quadratic term in b 



UJ 



tin 



+ - - m 2 + (m - b) 2 > 



(4.12) 



Notice that these are all positive, because £ > m. The minimum possible value, fixing £ and 
m but varying b, is given by b = m, in which case the frequency squared is 



to 



Im 



1 



m 



and the minimum value this can acquire is given by m = ±£, so that the mass is 

1 



^7 



+ 



(4.13) 



(4.14) 



So, as we go to higher and higher £, we find that the mode is more and more massive at 
the place where it is lightest (namely for b = m). However this only grows as which is 
subleading to the typical value of the frequency which is of order I. If we compare with our 



geometric result in equation (3.5) we find that it matches it very closely and it is exactly 



the same if we interpret (AL) 2 with the usual quantum value £(£ + 1) plus the 1/4 from the 
background curvature of the plane wave. 



4.2 Longitudinal fluctuations 

Now we analyze the fluctuations of the X 1 fields. This is trickier than for the Y a fluctuations. 
First, the X % have vevs, so that expanding the action in fluctuations is more involved. 
Secondly, the system is a gauged quantum mechanical system. This means that there are 
zero modes that should be projected out of the dynamics. Finally, for the displaced fuzzy 
spheres we are not at an extremum of the potential, so the gradient of the potential does not 
vanish. This means that the Hessian that determines quadratic fluctuations is not invariant 
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under non-linear field redefinitions, unlike in the case of an extremum of the potential where 
the Hessian is a symmetric tensor on the tangent space of the corresponding configuration 
point. This is potentially problematic for the removal of the zero modes, as the group of 
U(N) rotations of the configuration gives us a non-linear geometric space. 

All of these problems are solvable in practice. What we need to do is to argue that our 
fluctuations are orthogonal to linearized gauge transformations on a particular configuration 
defined by a background. Since we have a metric on the configuration space defined by the 
kinetic term, this is a well defined procedure. 

We expand the fluctuations of the off-diagonal blocks as follows^ 



* L3+b (° i)+E fa ^+(^Lr^L ) 

^ ' i,m 

x+ * L + + Y, 5x tm-iYem + (6xJ m+i yYl, 

X- * L +E 5 W^ + (^ + m -i)^!n- ( 4 -!5) 

Notice that we have shifted the index m in the coefficients of the X fluctuations by ±1. 
The reason for doing this is that the matrices L + and L~ usually are associated with angular 
momentum one. However the configurations with 5x = are spherically invariant, so the L + 
and L~ matrices should be associated to having no spin: in this way the spin of the matrix 
is cancelled by the spin of the label. The kinetic term for the longitudinal fluctuations will 
be then given by 

Notice that the metric on these fluctuations is diagonal, but the coefficients are not one. 
This is important for evaluating frequencies and for projecting out the gauge fluctuations. 

Now let us perform a gauge transformation that is off-diagonal, with parameters 59i m Yi m + 
h.c. This is necessary for the generator to be hermitian. We find that 

5 e X 3 = i[56 lm Y lm + h.c.,X 3 ] = -i(m - b)59 lm Y im + i(m - b)56* lm Y} m . (4.17) 

Again, notice that b just shifts the value of m for this calculation. Similarly we find that 

S e X + = -iy/(£ - m )(£ + m + l)6e lm Y lm+1 + iy/(e + m)(£-m + l)^ m ^l_i , 
5 e X~ = (S e X + y. (4.18) 

We then require that the allowed 5x\ m are orthogonal to the 59e m variations of the con- 
figuration. The conjugate variables to these rotations vanish. Classically this means that 



6 An equivalent way of doing this computation is to expand the fluctuations in the basis of eigenstates of 
the 6 = problem, which was originally solved in |21) . We outline this alternative derivation in the appendix. 
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for fluctuations proportional to S6e m we have to impose the constraint 58e m = 0, so that 

se tm = o. 

These gauge transformations are unphysical and are projected out by the Gauss law 
constraint. To match our labeling, we should replace the dummy index m — > m — 1 or 
m — > m + 1 in the various terms in the expansion of 5gX + and 5qX~ . We get this way that 

S e X + = -iy/(£ + m)(£- m + l)5e £m ^Y im + %y/{l -m)(£ + m+ l)5d} m+l Y^ m , 
5 e X~ = (5 e X + Y (4.19) 

and we see the consistency of the conventions of spin labeling, for the S9i m appear in the 
same way as the 5x\ m in the expansion of the fields. Notice how in the gauge variations the 
coefficient of Y& in 5qX~ vanishes, and also the one of Y^_^ in 5eX + . This is the object with 
maximum helicity (total spin along the 3-axis) at fixed t. It has helicity I + 1 and — i — 1 
respectively. From our previous considerations these are the modes that are most likely to 
become very light. We will see that they can indeed become tachyonic for some values of b. 

Now we move on to analyze the potential for these fluctuations. Using the following 
identities 

i[X\ X 3 ] + X 1 + % (t[X 3 , X 1 } + X 2 ) = [X+, X 3 ] + X+ , 
i[X 2 ,X 3 ]+ X 1 ~i (i[X 3 ,X 1 ]+ X 2 ) = -X~,X 3 ] + X~ , 

*[X\X 2 ] = i[X-,X+], (4.20) 
the potential can be rewritten as 

2 



V (X) - -tr 

"BMW _ tl 



1 



BMN 



2 



Q[x-,x+] +X 3 ) + ([X + ,X 3 ] + X + )(-[X~,X 3 ] +x- 



(4.21) 



If we expand in quadratic fluctuations, we can expand each term in the square to linearized 
order and get that way some quadratic terms. There is an additional term that arises because 
when we turn on b the background does no longer satisfy [X + ,X~] = 2X 3 . The other two 
equations are satisfied. These contributions from the potential not being at a minimum 
affect the coefficient of 5x + 5x~ in the quadratic terms. 



Writing (4.15) in blocks we have that, schematically, 



( ( .-V.T)< 

plus a similar expansion for 5X~ . To expand to quadratic order the total potential in off- 
diagonal modes there are two contributions: those that are linear in the fluctuations in 
each term of the potential that is squared and those that are quadratic in the fluctuations. 
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The linear terms are off-diagonal, the quadratic terms are block diagonal. The following 
intermediate calculations are useful before giving the final answer: 

5X 3 + 1 -[L-,5X + ]- 1 -[L + ,5X-] = 

Yl ( 6x *m + \^^-rn)(t + m + l)5x+ m - l -^{£ + m){£ - m + l)6xj m J Y im + h.c. , 



[5X + ,(X 3 )] + [L + ,5X 3 } + 5X + = 

J2((b-m + 1)5x1^ + y/(£ + m)(e-m + l)6xl m _^ Y tm 

+ ((m-b + l)(5xj m+1 )* - ^(£-m)(£ + m + l)(5x 3 tm+1 )*) Y} m , (4.23) 



and a similar equation involving 8X~ . All these terms are off-diagonal (being proportional 
to the Yi m ). Notice that because of our conventions, only the same values of £, m appear in 
all of the coefficients of these linear terms: this is, the mixing of modes only mixes the same 
values of £, m. This simplification makes the problem very tractable for these modes as well. 
In the end, we need to understand how three modes mix, but one such mode is projected out 
because of the gauge constraint. The general mass reduces to diagonalizing a 2 x 2 matrix 
for each value of £, m 

This can be combined to give a block diagonal term in ^[X~,X + ] + X 3 proportional to 

0\ 1 (8x-{8x-)*YY^ - 8x+(8x+)*YY* 
ftj + 2\ 5x + (5x + )*Y^Y -5x-(5x-)*Y^Y 

(4.24) 

When we square and take traces, expanding to quadratic order in fluctuations, this gives us 
a contribution to the mass matrix equal to 

\b{8x + {5x + Y - Sx-{Sx~)*) . (4.25) 

Notice that the contribution from this term is negative or positive for different modes de- 
pending on the sign of b. 

The kinetic term (4.16) suggest that we normalize the fields slightly differently, 8x = 
\/2 5X ± , to have canonical normalizations for every mode. The terms in the expansion can 
be rewritten as 



8x 3 m + -^(£-m)(£ + m+l)8x+ m - -^/(£ + m)(£ - m+ l)8xj m 

( Sx £m \ 

J (i-m)(e+m+ij _ ^ _^±^^£^ UXl m \ =V 3 5X . (4.26) 
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Similarly we find a V + and V given by 

V + = ^(£-m)(£ + m + l), 
V~ = (-y/(e + m)(£-m + l) 



V2(b - m) , 
0, V2(m 



b) 



(4.27) 



We have shifted m — > m ± 1 in the equations above so that we are comparing the same 
coefficients of £, m. 

The mass matrix is given by squaring these vectors and adding them together including 
also the contribution of (4.25) 



im 



(vyv 3 + ^{v + yv 



1 




(4.28) 



The end result is given by 

1 + I + I 2 - m 2 
= | (b-m+ 1)A_ 
(b-m - 1)A+ 



(6 - m + 1)A_ 
6 +(6 



(6 



A 2 



m - 1)A^ 
-A + A_ 



(4.29) 



-A,A_ 



-b +(b 



m) 



+ A 2 



where we have defined the shorthands 



A 



± 



(£±m)(£^m + l) 



(4.30) 



2 

Of particular interest to us is when m = £ + 1 for SX^ m (i.e. A + = 0), and when m = — £— 1 
for SXg m (i.e. A_ = 0). For these cases there is no mixing with any other mode and these 
fields have maximum spin in the 3-direction for fixed £. We have already argued why these 
modes are important. Their masses are given by 



-i-ij 



-b+(b-£-l) 2 
b+(b + £ + l) 2 . 



(4.31) 



These modes are tachyonic for b = ±(£+1) on an interval for b of order \f£. Notice that there 
is a tower of tachyonic modes for each b labeled by £ with a quadratic dispersion relation. 
This can be interpreted as a tower of tachyonic modes on a circle in the presence of some 
holonomy for a gauge field under which these fields are charged. Other modes for which m 
is not maximal in the sense above are not tachyonic. 

From the equation for the masses above we still need to project out the gauge variations. 
This is straightforward, but tedious. If we call the projection matrix that projects onto the 
gauge degrees of freedom as S^ m , then 1 — S£ m is the projection in the orthogonal components. 
The mass matrix we need is then given by 



Ira, phys 



Sfon)c^ m (l — S£ m ) 



(4.32) 
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The precise expressions are not very illuminating. However, none of the modes that appear 
this way are tachyonic, except the ones that we have already discussed. 

We can also check that the eigenvalues of the above matrix are £ 2 , (£ + l) 2 , when we 
set b = (as originally found in [21]) as a consistency check. For b = the modes with zero 
eigenvalue are the gauge zero modes. For b ^ these modes seem to become massive (the 
determinant is not zero), but as we have argued already, this is an artifact of the linearization. 
After all, expanding to second order in these gauge variations we find that 

vg} N (b, 56) ~ V™ N {b) + d h V^ N {b){89 2 ) + ..., (4.33) 

and the second term only vanishes for b = 0. However, the potential is invariant, so b must 
be corrected to second order in gauge fluctuations. This is a non-linear change of variables. 
This is why it is better to project on directions orthogonal to the gauge transformations than 
trying to sort this second order variation and how it affects the metric of the other modes. 

It is clear that when we consider the above result, we should organize the modes according 
to the following criteria. If the two fuzzy spheres intersect for some b > 0, we should fix 
K = £ — m for the modes where b ~ m (those that are near the intersection) . For each such 
K we get a tower of states on a circle labeled by the different values £ (or m) (there are two 
states for general £,m). These look like a tower of fields in one dimension. For b < we 
would do the same by fixing K = £ + m for the modes near the intersecting fuzzy spheres. 

Notice that this result is very similar to results that have been found in other matrix 
models with fuzzy spheres [33]. Since in their case the solutions with displaced fuzzy spheres 
are critical points of the potential, the issues with zero modes do not appear. 

4.3 Interpretation from J\f = 4 super Yang-Mills 

We would like to look at what we have obtained so far from the point of view of M = 4 
super Yang-Mills on S 3 x R. The results above are expected to give details of an SU{2) L 
invariant subsector of this theory on the sphere. After all, the BMN dynamics is exactly 
the truncation to the SU(2)l invariant subsector of the field theory. Thus, any dynamical 
feature present in the SU(2)l invariant truncation is a feature of the full M = 4 super 
Yang-Mills on the sphere. In particular, the unstable modes we found in the matrix model 
are necessarily unstable modes in the full field theory. 

For all the fluctuations we have studied, notice that, apart from the bounds on the total 
magnetic moment, \ji ± j 2 1 , our results are essentially independent of the values rii and n 2 . 
If we think back at how the plane wave matrix model is related to a truncation of M = 4 
super Yang-Mills on S 3 we realize that we should not be surprised by this fact. The splitting 
into n\ and n 2 as fuzzy spheres is artificial in M = 4 super Yang-Mills. The fuzzy spheres 
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are gauge transformations of the vacuum. However b matters: its presence gives states with 
non-zero energy. In the field theory setup the total angular momentum of physical states also 
matters, but not the details of the splitting into n 1; n 2 , as these would just give a multiplicity 
of components with different masses from different angular momentum objects. 

To make things precise, we should remember that the fields X % arise from the gauge 
connection on S 3 , while the fields Y a arise from the constant scalar modes of the field 
theory. The way to do this is simple: since we have an S 3 sphere, the spatial components of 
the gauge connection can be written in a basis of left invariant one-forms e l 

A(t) =X l (x,t)e i . (4.34) 

The spatial component of the field strength (chromo-magnetic field) is given by 

dA + A A A = dX l Ae' + XW + - [X\ X^e* A e j . (4.35) 

Restricting to SU(2) invariant states requires that X 1 be constant. Thus on such configu- 
rations the first term vanishes. However, the one-forms are not exact de % 7^ 0. Instead they 
satisfy the Maurer-Cartan equations. Thus the second and third terms do not vanish. Even 
for constant diagonal X 1 there is a non-zero result for the magnetic field. Considering the 
magnetic field squared part of the Hamiltonian gives an expression that is exactly the term 



in the potential given by equation (2.6), after one accounts for various normalization issues. 
Thus, a non-trivial diagonal displacement in the X % (of size b) corresponds to a chromo- 
magnetic field of strength proportional to b (the commutator squared term vanishes for this 
configuration) . 

Notice that in such a description of the theory the parameter b would correspond to 
a magnetic field that leaves invariant an SU(2)l of rotations on the S 3 (the electric field 
is related to the time derivative of the X 1 ). It also leaves invariant an U(1)r of rotations 
along the direction of said magnetic field. For such a system the SU(2) spherical harmon- 
ics under SU{2) L and the £7(1) quantum numbers completely specify the spectrum for the 
scalar fluctuations, after all, these are originally classified as the (n, n) representations un- 
der SU(2) x SU(2). The first quantum number determines n, while the second quantum 
number denotes the ^-component of spin in the SU(2)r. Moreover, this symmetry splits the 
degeneracy in energies for fluctuations with different U(1)r completely. 

When truncating to SU(2)l invariant states the choices of fuzzy spheres matter again, 
but this does not affect the energies, just the SU{2) L labels of the fields that get twisted. 
As we also argued, the SU{2) R labels are tied very closely to the SU{2) L labels of the 
field around the trivial vacuum. Since the system preserves the SU(2)l symmetry, we find 
that the splitting of modes into SU(2)r representations is preserved in the trivial vacuum. 
This is not obvious in the matrix model computation since we were turning on a vacuum 
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expectation value that breaks the SU(2)r symmetry to U(1)r and in principle could induce 
mixing between modes with different SU(2)r representations. This non-mixing between 
modes looks like a happy coincidence from the matrix quantum mechanics point of view. 
Here we see that there is a symmetry reason from super Yang-Mills for that to happen. 
Moreover, each such (n,n) can contribute at most one singlet under SU(2)l for each value 
of the £7(1) angular momentum. 

We should also ask what kind of instability do the modes between intersecting fuzzy 
spheres represent in super Yang-Mills. Clearly, the fuzzy sphere is a gauge artifact, but not 
the presence of a magnetic field. So we should ask what kind of instabilities can arise in the 
presence of a constant chromo-magnetic field on a sphere. If we take b large, this corresponds 
to a large magnetic field. By thinking in terms of engineering units, the magnetic field b 
generates a scale in the system that, if b is large, is much shorter in length than the radius 
of the sphere. Under such conditions we can ignore the radius of the sphere and treat the 
magnetic field as if it were constant in space. The only scale in the problem is associated 
to b. The mass squared of the unstable mode should therefore be linear in b (just a result 
of dimensional analysis). If one has momentum p along the direction of the magnetic field, 
the spectrum gives modes with frequencies of the form ujf m = —\b\ + p 2 . This is exactly the 
result we found in equation (4.31 ), if we identify p = £ + 1 ± b. 

Such a result is well known. It is the Nielsen-Olesen instability [31] (see also [35]). This 
instability is due to the fact that gluons charged under the field that acquired a magnetic 
field have a large magnetic moment that is enough to drive them tachyonic. For scalars there 
is no magnetic moment, and the localization effects in a magnetic field cost energy via the 
uncertainty principle, so these modes are not tachyonic, just as we found for the fluctuations 
of the Y a fields. For fermions, one can get zero modes, which is the familiar connection 
between massless fermions and index theory. Incidentally, it is well known that it is the 
magnetic moment contribution to the /3-function of non-abelian gauge theories that turns 
the sign contribution for vector particles relative to scalar particles. It was argued that this 
would also lead to an effective action where the chromo-magnetic field condenses j3S]- The 
Nielsen-Olesen instability would destroy that type of order. 



5 Time dependence for the tachyonic modes 

So far in our analysis we have considered configurations that are frozen in time. We have 
seen that tachyons can form in the region where the two fuzzy spheres intersect and we have 
argued that they do not mix with any other mode. At this point an interesting question to 
ask is what happens to these tachyons as we let our system evolve. In this section we present 
a straightforward analysis for the simplest possible trajectory, namely an oscillation of the 
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spheres along one fixed direction. 

Consider two fuzzy spheres of sizes ji and j%. The tachyons appear for any integer value 
from \ji — J2I all the way to j\ + j'2 if the smaller sphere is allowed to oscillate from the center 
all the way to the outside of the larger sphere. As we have computed in (4.31), the mass for 
these tachyons is given by 

(mf(t)) 2 = (£+l±b(t)) 2 ±b(t), (5.1) 

where we do not write the label m, which is set equal to =f(£ + 1). We take b(t) = frsint, so 
to have a simple sinusoidal motion along the X 3 direction. The motion in t is periodic with 
period 2n, so the full analysis can be restricted to the interval t G {0, 2n}. We need to solve 
the differential equation 

q e (t) + (mf(t)) 2 q(t) = 0. (5.2) 

This is the equation for a harmonic oscillator with time dependent mass. It is somewhat 
similar^] to the Mathieu equation that describes parametric resonances and, in cosmology, 
the fluctuations of the inflaton around the minimum of the potential during preheating [18J. 



In general there are two linearly independent solutions to the equation (5.2), which we 
call qi(t) and (fe(£)- We can relate the initial time problem (at time t) to the problem at one 
period later (at time t + 2ir) using a periodicity matrix 



qi (t + 2n)\ = (A B\ ( qi {t) 
q 2 (t + 27r)J \C D) \q 2 (t) 



(5.3) 



This equation can be diagonalized, so we can choose the solutions to be eigenvalues of the 
matrix above. The Wronskian of the solution is constant, so the matrix transforming between 
one and the other has determinant equal to one. Also, since the differential equation has 
real coefficients, the solutions can be made real and in that case the matrix above is real as 
well. Hence, the eigenvalues are either real or unitary. These eigenvalues serve as Lyapunov 
exponents for the classical periodic orbit. When the eigenvalues are unitary the system is 
stable, when the eigenvalues are real the system is unstable. 

The system can be interpreted also as a Schrodinger problem with fixed energy in a 
periodic potential]^] which is the negative of the (mf(t)) 2 function. If the solutions are 
quasi-periodic (the eigenvalues are in the unit circle), one of the functions is identified with 
positive frequency and the other one is identified with negative frequency modes. This is 
the case where the functions qi(t) belong to a band of the periodic potential. 

Generically, if there are regions where the mode is tachyonic, the corresponding Schrodinger 
particle needs to tunnel through the barrier. This phenomenon generically leads to the prop- 
erty that the eigenvalues of the matrix above are non-unitary, with the tunneling amplitude 



7 Albeit more complicated and not generically solvable in terms of elementary functions. 
8 For an elementary treatment see [37] . 
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characterizing the growth of the signal. We can in general estimate this using a WKB 
approximation. 

The large eigenvalue tells us how the modes grow around these periodic solutions and it 
describes the discrete time dependence of the instability under various oscillations. The two 
linearly independent solutions can also be thought of as coefficients of raising/lowering oper- 
ators. The matrix computed in this basis is a Bogolubov transformation for each period and 
the amplitude growth correlates with the amount of particle creation between oscillations. 

For us, the most important question to ask is which of all the modes above grows the 
fastest, as these modes will dominate the initial stages of the brane collapse problem. It is 
not hard to figure out that the tachyon with the highest i will generally dominate. First, it 
will be tachyonic for the longest time during the periodic trajectory, not only because the 
range of b where it is tachyonic is larger (recall that this range is of order \fl), but also 
because the motion of the oscillations is slower at larger displacement. This means that the 
main condensation of modes will happen between the north pole of one fuzzy sphere and the 
south pole of the other one. We can do this numerically for various values of the amplitude 
of oscillation b and for different values of £. This is shown in Figures [3] and H 




Figure 3: Here we show the norm of the maximum eigenvalue of the periodicity matrix for 
various values of I as a function of the amplitude b, labeling the horizontal axis. We clearly 
see the band structure of the problem and that the largest value of i is typically the one 
with most amplification so long as it is tachyonic. 
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Figure 4: Here we show the maximum eigenvalue of the periodicity matrix for various 
values of £, where we explore the problem in the region of large amplitude b, labeling again 
the horizontal axis. Asymptotically the amplification becomes of order one, since the time 
during which each mode is tachyonic on a single oscillation becomes very small. 



Notice also that because condensation happens for various modes with different £, the 
system classically breaks the rotation symmetry around the axis of symmetry for small 
perturbations on the off-diagonal modes. Each of these is expected to decohere from the 
others once the non-linearities set in, because the classical trajectories diverge rapidly from 
each other. Also, as seen in the figures, any small amount of back-reaction can move modes 
between the oscillation bands and the growth bands. Let us assume that this will happen 
for some fixed value of the off-diagonal perturbation in the classical picture. The time to 
reach this value will depend on the initial amplitude of the off-diagonal mode. If we begin 
in an adiabatic ground state for these modes in the non-tachyonic region, we can estimate 
the initial size of fluctuations by a harmonic oscillator wave function. This is a power-law 
function of h. The growth is exponential (and can easily be of a few orders of magnitude per 
oscillation), so the time growth for quantum perturbations around a fixed classical periodic 
orbit to the stages where back-reaction becomes important is logarithmic in the ft — )• limit. 
As we increase h, the time to reach back-reaction is shorter. In the very quantum regime 
the approximations used here break down and a more robust formalism with quantum back- 
reaction needs to be developed (an initial attempt of such a setup has been recently pursued 
in [2B])- This suggests that in the strong coupling regime thermalization could be very fast. 
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6 Conclusion 



The purpose of this paper was to initiate a study of the process of formation and thermal- 
ization of black holes, using the gauge/gravity correspondence. Since these are notoriously 
complicated problems to address in ordinary field theories, our strategy has been to focus 
our attention on simplified models with a reduced number of degrees of freedom. 

The prototypical example of such models is the BFSS matrix quantum mechanics de- 
scribing M-theory on the discrete light-cone quantization of flat space. Unfortunately, this 
model possesses some features that make its study extremely challenging. These features 
include the presence of flat directions in the potential (which give rise to a continuous spec- 
trum and to the difficulty of distinguishing between single-particle and multi-particle states) 
and the absence of a tunable coupling constant. Moreover, the wave function of bound 
states in this model is not known, making it impossible to describe scattering processes in 
the regime of high energy and small impact parameters where the non-linearities of gravity 
become important and black hole formation takes placer^ 



We have therefore considered a different model, the BMN matrix model given by (2.4), 
where many of the problems that plague the BFSS model no longer occur. For example, the 
mass terms in the BMN model lift flat directions and give rise to a discrete spectrum with an 
isolated set of classical vacua, the fuzzy spheres that have been the central ingredient of our 
setup. In the limit of large mass, these different vacua are divided into superselection sectors 
described by harmonic oscillators and the spectrum of their fluctuations can be studied. 
While the supersymmetric vacua of the model consists of concentric spheres, one can also 
consider non-BPS configurations obtained by displacing the centers of the spheres. These 
are obtained by turning on the diagonal modes that control the center of mass motion of 
the spheres without deforming their shape. In this way one can setup a scattering problem, 
that might be used as a proxy, under computational control, for the high energy scattering 
of particles in gravity. 

We have started our analysis by considering two fuzzy spheres, displaced along a direction 
in their world-volume and frozen in time in that position. The formation and thermalization 
of a black hole can be associated to copious particle production in off-diagonal modes of 
a configuration. The presence of classical tachyons makes the analysis simpler as we get a 
classical instability that can in principle drive the system towards thermalization. Expand- 
ing the fields in fuzzy spherical harmonics, we have found that the modes with the maximal 
angular momentum along the direction of the displacement become tachyonic at the inter- 
section locus between the two spheres. Interestingly, this instability has a four-dimensional 

9 Most of the work in the past aimed at matching the loop expansion of the BFSS matrix model with 
perturbative expansions in gravity has been mostly limited to linearized examples or to terms protected by 
non-renormalization theorems, see for example |39j and (40] for a review. 
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interpretation. It can be regarded in fact as a Nielsen-Olesen instability in M = 4 super 
Yang-Mills, of which the BMN model is a truncation. The role of the background magnetic 
field is played in our system by the displacement vector. 

An obvious generalization of our initial conditions would be to allow for a displacement 
of the fuzzy spheres also along the transverse directions (this corresponds to turning on 
some (Y a ) ^ 0). The analysis in this case is slightly complicated by the fact that the 
fluctuations of the X 1 and the Y a fields get coupled and the diagonalization problem becomes 
less straightforward. One could also try to include the fermionic modes, which have not been 
considered here. 

In the last part of the paper, we have allowed our initial configuration to evolve in time 
describing a periodic motion, with the two fuzzy spheres oscillating along the direction of 
the displacement and crossing each other repeatedly. This is the simplest trajectory to 
study and gives rise to a dynamics that is somewhat similar to the physics of preheating 
during inflation. We have argued that the tachyonic modes that form at the intersection 
locus between the spheres typically get reinforced after each period of the oscillation. The 
growth is exponential, and can potentially give rise to a fast thermalization of the other 
degrees of freedom living on the spheres. So if these instabilities would cause the BMN 
model to thermalize fast, it is suggestive that the Nielsen-Olesen instability could drive fast 
thermalization in QCD processes like heavy-ion collisions. This thermalization is observed 
in experiments at RHIC [H] and it is argued to happen in a very short time scale [42J. 

A better understanding of the details of this dynamics is surely desirable and this paper 
can be considered as a first step toward this ambitious goal. In particular, we find it extremely 
interesting to try to estimate from our setup the time scales characterizing the various phases 
of the black hole evolution, for example the thermalization time. To this regard, it has 
recently been conjectured [I] (prompted by the analysis in [13]) that black holes are the 
'fastest scramblers' in nature. By fastest is intended that the thermalization time scale is 
logarithmic in the number of degrees of freedom of the system, rather than a power-law, as 
was originally proposed in |H]. Our hope is that it might be possible to check this claim 
using the ideas and techniques we have presented in this paper. 

An even more involved scenario would start with a configuration that carries angular 
momentum on the plane. Such a scenario might lead to different black hole shapes. The 
time dependent analysis we performed would be more complicated because there is more 
mixing between modes (the system does not preserve an azimuthal symmetry). 

We conclude by repeating the observation in the Introduction that the BMN matrix 
model is amenable to being put on a computer. We can think of implementing numerical 
simulations of our system, where we define some initial configuration of displaced fuzzy 
spheres and let them evolve. Such simulations should provide us with a more accurate 
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description of the details of the dynamics that follows the formation of the tachyons and 
shed light on what happens during the first phases of the evolution of the black holes. We 
are currently looking into this. 



Acknowledgements 

D.B. would like to thank C. Asplund, D. Kabat, J. Maldacena, E. Silverstein, and H. Verlinde 
for various discussions related to this work. D.T. is grateful to S. Giddings for several 
instructive discussions on high energy scattering in gravity and to N. Iizuka for collaboration 
on related topics. D.B. would like to thank the Simons Center for Geometry and Physics, 
where some of this work was carried out. Work supported in part by DE-FG02-91ER40618. 
D.T. was also supported by PHY04-56556 and DE-FG02-95ER40896. 



A An alternative derivation of the longitudinal spec- 
trum 



In this appendix we present an alternative derivation of the spectrum of longitudinal fluctu- 



ations (4.31). The idea is to expand the fluctuations in the basis of eigenstates of the 6 = 
system. 



We start by rewriting the quadratic potential for the 5X l fluctuations as 



V, 



1 



mtx - -tr ( {6X* + i^ k [(Xt),6X k ]) -iWX* 




b 



5X 1 



(A.l) 



The case with 6 = was worked out in detail in [21], where it was found that the solution 
to the eigenvalue problem 



SXL.+iP" [L j ,5XU =X a 5X t 



(A.2) 



is X a = (0, -£,£ + 1) with 



(£-m)(£ + m + l) ^ 

5X (o) = JuTT) +1 + 



(e + m)(e-m + l) 

0A -(0) ~ r Yim\l mn , 1N 1 £m-l + n - c - > 



£(£+1) 



m 



Ytm + h.c. 
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SX, 



i-t) 



sx, 



(-£) 



sx. 



SX, 



(M-l) 



— Oi£-l,r 



at-i. 



(£ + m)(£ + m + V 



£(2£ + 1) 



Y lm+ x + h.c. 



'(£-m)(£-m + l) 



£(2£ + 1) 



Ye m -i + h.c. 



\£ + m)(£-m) 



Y lm + h.c. 



(£+1) ~ A+i 



\£-m)(£-m + 1) 
{£+l){2£+\) 

\£ + m)(£ + m + I) 
{£+l){2£+l) 



Y lm+ i + h.c. 



Ye m -i + h.c. 



'(£ + m + l)(£-m + l) 



Y im + h.c. 



(£ + l)(2£+l) 

Expanding the off-diagonal ffuctuations in the basis of eigenstates above 

5X* = J2SX{ a) , 



(A.3) 



(A.4) 



the potential (A.l) becomes 



le 



»j3 




b 



5X U 




b 



We can write the blocks explicitly 

\tr £ (a q a, (our+ )t*x+ + (SX^SX^ + 2(SXf a) )UX i 



sx (0) 



(A.5) 



+b (b + \ a + \p + 1) (5X+ ) ) t 5X+ ) + b (b - X a - X p - 1) (5X M y5X i 



(ay U ^W) 



(A.6) 



Taking the trace, it is easy to see that the first line of the expression above reproduces the 
6 = spectrum 

\Y,P\<xi-i, m \ 2 + {l + 1) 2 |/WI 2 , (A.7) 



with the ■y£ m modes being the gauge variation zero modes. The second line gives instead 
(here we call 8xi a \ whatever is multiplying the spherical harmonics in (A.3), including the 
square roots) 



a,/3 £,m 
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+ (b - X a - X p - 1) (Sx (a) )* im (5x m ) tm 



(A. 



Expanding the sums in a and f3 in (A.8) and combining this with (A.7), it is straightforward 



to check that the resulting mass matrix admits the same eigenvalues obtained in (4.31). 
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